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Theory  and  Experimental  Validation 

A.  Bellettini,  M.A.  Pinto 


Executive  Summary: 

Higher  resolution  sonars  will  be  required  to  detect  and  more  importantly  to 
classify  mines,  which  are  undetectable  by  existing  sonars,  for  reasons  of  shape, 
material,  size  or  location.  Synthetic  aperture  sonar  (SAS)  will  increase  the 
sonar  cross-range  resolution  by  several  orders  of  magnitude  while  maintaining 
or  increasing  the  area  search  rate  and  thus  contribute  to  an  enhancement  of 
mine  hunting  performance,  in  particular  in  shallow  water,  where  smaller  mines 
are  more  effective.  The  side  scan  configuration  of  SAS  makes  it  well  suited  to 
remotely  operated  underwater  vehicles,  such  as  those  forseen  for  covert  survey 
and  reconnaissance  operations. 

SAS  performance  is  limited  by  the  precision  with  which  the  motion  errors  of  the 
platform  can  be  estimated.  The  terminology  of  “micronavigation”  is  used  to 
describe  this  very  specific  requirement  for  sub-wavelength  short-term  relative 
positioning.  The  aim  of  this  work  is  to  quantify  the  theoretical  performance 
of  data-driven  micronavigation  based  on  the  displaced  phase  centre  antenna 
(DPCA)  concept,  extending  the  results  obtained  in  Saclantcen  memorandum 
SM-352.  The  unique  feature  of  DPCA,  of  great  practical  significance,  is  that 
it  does  not  require  the  presence  of  seafloor  features  as  it  exploits  the  spatial 
coherence  properties  of  seafloor  reverberation. 
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Abstract:  The  Cramer-Rao  Lower  Bounds  (CRLBs)  on  the  cross-track 

translation  and  rotation  of  a  Displaced  Phase  Centre  Antenna  (DPCA)  in  the 
slant  range  plane  between  two  successive  pings  (known  as  DPCA  sway  and  yaw 
in  what  follows)  are  computed,  assuming  statistically  homogeneous  backscatter. 
These  bounds  are  validated  using  experimental  data  from  a  118-182  kHz  sonar, 
showing  an  accuracy  of  the  order  of  20  microns  on  the  ping-to-ping  cross-track 
displacements. 

Next,  the  accuracy  required  on  the  DPCA  sway  and  yaw  in  order  to  achieve 
a  given  SAS  beampattern  specification,  specified  by  the  expected  SAS  array 
gain,  is  computed  as  a  function  of  the  number  P  of  pings  in  the  SAS.  Higher 
accuracy  is  required  when  P  increases  to  counter  the  accumulation  of  errors 
during  the  integration  of  the  elementary  ping-to-ping  estimates:  the  standard 
deviation  must  decrease  as  P_1// 2  for  the  DPCA  sway  and  P~3/2  for  the  yaw. 

Finally,  by  combining  the  above  results,  the  lower  bounds  on  DPCA  micronav¬ 
igation  accuracy  are  established.  These  bounds  set  an  upper  limit  to  the  SAS 
length  achievable  in  practice.  The  maximum  gain  Q  in  cross-range  resolution 
achievable  by  a  DPCA  micronavigated  SAS  is  computed  as  a  function  of  the 
key  SAS  parameters.  It  is  found  that,  for  P  1,  the  optimum  SAS  spatial 
sampling  factor  is  4,  in  the  sense  that  it  allows  maximum  Q.  These  theoretical 
predictions  are  compared  with  simulations  and  experimental  results. 
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1 

Introduction 


The  performance  of  Synthetic  Aperture  Sonar  (SAS)  on  ocean  going  platforms  is 
limited  first  and  foremost  by  the  accuracy  with  which  the  motion  of  the  platform 
can  be  estimated  [1].  The  term  “micronavigation”  is  used  here  to  describe  this 
very  specific  requirement  for  sub- wavelength  short  term  relative  positioning,  which 
is  beyond  the  scope  of  instrumentation  for  most  high  resolution  imaging  applica¬ 
tions.  In  recent  years,  data  driven  micronavigation  techniques  have  emerged  as  a 
possible  solution.  The  most  promising  are  based  on  the  concept  of  Displaced  Phase 
Centre  Antenna  (DPCA)  which  exploits,  in  a  unique  way,  the  spatial  and  temporal 
coherence  properties  of  the  seafloor  backscatter  [2].  DPCA  is  a  known  concept  in 
radar  space-time  processing  [3].  DPCA  also  forms  the  basis  of  correlation  sonar, 
a  subject  which  has  recently  been  revisited  theoretically  by  Doisy  [4]  who  derived 
the  accuracy  of  translational  displacement  estimates  for  volumetric  arrays,  as  well 
as  for  attitude-stabilized  planar  arrays.  Several  experimental  SAS  prototypes  ex¬ 
ploit  the  DPCA  principle  for  SAS  micronavigation  [5,  6,  7,  8].  A  large  data  set  has 
been  collected  with  a  100  kHz  towed  SAS  to  evaluate  DPCA  micronavigation  for  a 
wide  range  of  parameters  [9].  The  aim  of  this  work  is  to  determine  the  theoretical 
accuracy  of  DPCA-based  SAS  micronavigation. 

In  sections  2  and  3  the  basic  principles  of  DPCA  are  summarized.  Next,  in  section 
4,  the  CRLBs  of  the  DPCA  sway  and  yaw  are  computed  and  compared  with  exper¬ 
imental  results.  In  Section  5,  the  accuracy  required  on  these  quantities  in  order  to 
achieve  a  given  expected  SAS  array  gain,  is  computed  as  a  function  of  the  number 
P  of  pings  in  the  SAS.  In  Section  6  the  results  of  Sections  4  and  5  are  combined. 
An  effective  reverberation  to  noise  ratio  peff  is  introduced  and  the  maximum  gain  Q 
in  cross-range  resolution  achievable  by  a  DPCA  micronavigated  SAS  is  computed 
as  a  function  of  peff  and  a,  the  SAS  spatial  sampling  factor.  It  is  also  shown  that, 
for  P  >  1,  a  =  4  is  optimum,  in  the  sense  that  it  allows  maximum  Q  for  given 
peff.  These  theoretical  predictions  are  compared  with  simulations  and  experimental 
results. 
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2 

Phase  Centre  Approximation 


Let  T  be  the  position  of  the  transmitter  centre  at  a  given  sonar  ping  and  R  the 
position  of  the  receiver  at  the  instant  of  the  reception  of  the  echo  from  an  arbitrary 
scatterer  located  at  X  (Fig.  1).  The  basis  of  the  Phase  Centre  Approximation  (PCA) 
is  to  replace  the  true  bistatic  situation  by  a  fictitious  monostatic  one,  which  assumes 
transmission  and  reception  occurs  from  C  —  (T  +  R)/ 2. 

To  see  when  this  is  valid,  we  compute  the  excess  in  round  trip  travel  path  for  the 
bistatic  case  over  the  monostatic  one.  Let  r  =  CX,  A  =  TR.  Expanding  in  series  of 
A/r  <C  1,  one  obtains: 

a  2  A  4 

TX  +  XR.  -  2  CX  =  —  cos2  6  +  — 7T  cos2  0(4  -  5  cos2  9)  +  . . .  (1) 

4  r  64 r6 

whore  9  is  the  bearing  of  X  (Fig.  1). 


X 


Figure  1  Geometry  of  the  Phase  Centre  Approximation.  T,  R,  C  and  X  are  re¬ 
spectively  the  positions  of  the  transmitter,  receiver,  phase  centre  and  single  scatterer. 

It  is  seen  that  PCA  holds  when  A2/4r  <  A0,  or  equivalently  r  >  A2/4Ao,  which  can 
be  interpreted  as  a  far  field  condition.  More  generally,  PCA  also  holds  in  near  field  at 
the  condition  that  the  received  signal  is  advanced  by  A2/ (4rc)  and  the  transmission 
sector  satisfies 

■X- (1  -  cos20e)  <  1  (2) 

4rA0 

where  9„  is  the  half  transmission  beamwidth.  This  condition  ensures  that  the  excess 
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round  trip  travel  path  is  the  same  for  all  scatterers  at  range  r  within  the  3  dB 
transmission  beamwidth,  so  that  their  interference  pattern  is  also  the  same. 

In  addition,  this  derivation  has  assumed  that  the  bistatic  angle  A /r  is  small  com¬ 
pared  to  the  beamwidths  of  both  the  transmitter  and  the  receiver,  so  that  the 
corresponding  changes  in  transmission  and  reception  directivity  gains  between  the 
bistatic  case  and  the  monostatic  one  are  negligible.  This  condition  is  almost  always 
valid  for  usual  SAS  systems. 


-3- 


SaCLANTCEN  SR-355 


3 

Displaced  Phase  Centre  Antenna 


The  following  analysis  of  DPCA  applies  to  a  multi-element  SAS,  which  consists  of 
a  transmitter  of  length  L(  =  d  and  of  a  receive  aperture  of  N  >  1  elements  spaced 
at  d,  of  total  length  Lr  =  Nd.  Let  Tp  be  the  position  of  the  transmitter  centre  at 
ping  p  and  RnP  that  of  receiving  element  n  at  the  instant  of  the  reception  of  the 
echo  from  range  r.  Based  on  PCA  this  is  equivalent  to  transmission  and  reception 
from  Cnp  =  (Tp  +  Rnp)/2. 


The  basic  idea  of  DPCA  is  to  cancel  the  along-track  displacement  of  the  sonar 
between  two  successive  pings  by  synthesis  of  an  effective  displacement,  in  the  op¬ 
posite  direction,  of  a  subset  of  receiving  elements.  This  is  achieved  by  operating  at 
D  =  Md/2  (. M  integer),  with  M  <  N  so  that,  in  the  absence  of  cross-track  motion, 
there  ar e  K  =  N  -  M  phase  centres  which  overlap  from  ping  to  ping  (see  Fig.  2) 
and  which  form  the  DPCA.  It  follows  that  the  signals  of  the  DPCA  are  identical 
for  both  pings,  up  to  noise  and  possible  changes  in  the  propagation  medium  or  the 
scatterer  geometry. 

In  the  presence  of  cross-track  motion,  the  cross-track  displacements  7 k  of  the  K 
DPCA  elements  are  given,  for  a  small  angle  ip,  by 

Ik  =  7  +  V><4  (3) 

where  7  and  ip  are  defined  as  the  DPCA  sway1  and  yaw  (Fig.  3),  and 

is  the  abscissa  of  the  phase  centre  of  DPCA  element  k,  with  respect  to  the  origin 
taken  at  the  centre  of  the  DPCA. 

The  cross-track  motion  leads  to  a  change  in  round-trip  travel  path  to  scatterer  X 
equal  to  27^.  cos  0,  where  6  is  the  bearing  of  X.  Under  the  condition  that 

^(l-cos0e)<l  (5) 

Ao 

LIt  should  be  cautioned  that  the  DPCA  sway  differs  from  that  of  the  physical  array  centre  in 
the  presence  of  yaw.  The  consequences  of  this  will  be  explicated  in  Section  5. 
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the  round-trip  travel  path  to  all  the  scatterers  within  the  transmission  beam  is  the 
same  and  equal  to  2 7*,.  It  follows  that  the  DPCA  signals  will  again  be  identical  for 
both  pings,  up  to  a  delay 

Tk  =  27  k/c  (6) 

which  can  be  estimated  by  a  cross-correlation  over  a  short  range  window  centred  at 
r. 


L  =  Nd 

_  _  J 

[ .  -TJ 

(N-M)d 

Rip*  1  •  lRup  • 

1  *  |FV* 

d/2 


N-M  overlapping  phase  centres 

CN0Ht) 


r~^ 


I  *  |R(n-»Vii 

• 

!RN(fH-lT  1 

<D  =  Md/2> 

,  (N-M)d  t 

Ping  p+1 


Figure  2  Displaced  Phase  Centre  Antenna  in  the  absence  of  cross-track  motion 
(the  upper  array  has  been  shifted  for  clarity). 


Figure  3  DPCA  sway  7  and  yaw  4>  between  successive  pings. 


It  should  be  noted  that  condition  (5)  limits  the  maximum  cross-track  motion.  When 
(5)  no  longer  holds,  the  DPCA  signals  decorrelate.  This  effect  has  been  studied  in 
detail  for  interferometric  sonar,  where  it  is  known  as  “baseline  decorrelation”  [10]. 

Furthermore,  it  has  been  assumed  above  that  the  along-track  sonar  displacement 
was  known  and  equal  to  Md/2  with  M  integer.  A  standard  navigation  system  could 
be  used  for  this  purpose,  together  with  a  system  that  modifies  the  ping  repetition 
period  to  attain  the  desired  spatial  sampling,  since  the  accuracy  requirements  are 
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not  as  great  as  for  cross-track  motions.  Alternatively  DPCA  can  be  extended  to 
estimate  these  quantities  as  well.  However  this  will  not  be  discussed  further  here. 

By  integrating  the  ping-to-ping  displacement  estimates  provided  by  the  DPCA,  the 
projection  of  the  platform  trajectory  in  the  slant-range  plane  can  be  reconstructed 
and  used  to  form  the  SAS  image.  Clearly,  the  finite  accuracy  of  the  DPCA  will  limit 
the  number  of  pings  P  which  can  be  summed  coherently,  possibly  to  a  fraction  of 
the  maximum  number  determined  by  the  transmission  beamwidth  and  range.  The 
corresponding  cross-range  resolution  gain  Q  achievable  in  this  manner  is 

rr  +  2(r-i.).g  =  1  +  pzi  (7) 

hr  a 

where  a  =  Lr/2D  is  the  SAS  spatial  sampling  factor.  The  problem  posed  is  the  eval¬ 
uation  of  the  accumulated  error  of  DPCA  micronavigation  and  the  limits  imposed 
on  the  achievable  Q. 
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4 

Cramer-Rao  lower  bounds  on 
DPCA  sway  and  yaw 


4.1  Theory 


Let  {Xi(t),  X[(t), X2(t),  X^it),  ...,Xx(t),X'K(t)}  be  the  K  pairs  of  signals  received 
by  the  elements  of  the  DPCA  at  two  successive  pings.  We  will  assume  that,  over 
the  temporal  estimation  window  of  duration  W,  these  signals  can  be  expressed  as 

Xk(t)  =  Sk(t)  +  Nk(t)  (8) 

X'k{t)  =  Sk(t-Tk)  +  N'k(t)  (9) 

where  the  noises  Nk  and  Nk  are  mutually  independent  and  independent  of  the  re¬ 
verberation  signals  Sk.  Furthermore,  the  Sk  are  assumed  independent  from  each 
other.  This  holds  because  adjacent  elements  are  separated  by  d  =  X(,  which  is  the 
spatial  correlation  length,  for  statistically  homogeneous  reverberation  [4].  Finally, 
reverberation  and  noise  will  be  assumed  to  be  Gaussian  random  processes  with  a 
flat  power  spectral  density.  The  reverberation  to  noise  ratio  will  be  denoted  by 
p  =  (S*)/(N2). 

To  compute  the  CRLBs  we  proceed  to  compute  the  log-likelihood  relative  to  the  two 
parameters  (7,^)  and  the  corresponding  2x2  Fisher  Information  Matrix  (FIM).  As 
(A/c,  A£)  and  (A/,  X[)  are  independent  for  k  ^  l,  this  log-likelihood  can  be  expressed 
as  a  function  of  log-likelihood  C\ (D)  relative  to  the  estimation  of  a  single  delay  D 
as  follows: 

£k{ 7,V0=  5Z  A(rfc(7,V’)).  (10) 

1  <k<K 

The  required  FIM  F(7 ^  can  now  be  obtained  directly  as  a  function  of  K  and  the 
FIM  of  order  1  relative  to  the  estimation  of  a  single  delay  t,  defined  as 

Ft  =  -  (d2A/dr2) .  (11) 


One  has 


±/££l\(  k  E dk\ 

c1  \  dr2  j  \Yjdk  E  4  ) 
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4  (  K  0  \  /10, 

°  0  )  ■  (12) 

The  zero  non-diagonal  terms  show  that  the  estimations  of  y  and  %jj  are  separable. 
The  expression  of  aT  =  is  a  result  of  time  delay  estimation  for  passive  sonar 

(see  (23)  of  [11]):  _ 

CTt=  27t/oV/w/p  + 2p2  ^ 

where  B  is  the  signal  bandwidth. 


The  CRLBs  can,  after  some  elementary  manipulations,  be  written  as: 


1  A0  1  1  /II 

~  2tt  2  sfBW  s/K  V  P  V 

and 

y/Z  Ik -7  Ao  1  1  /l  1 

+  1  (K-l)d  JBW  v^V  P  v* 


(14) 


(15) 


The  physical  interpretation  of  these  formulas  is  straightforward.  The  CRLB  on  y  is 
proportional  to  Ao/2,  the  half- wavelength  at  the  centre  frequency.  The  CRLB  on  ^ 
is  proportional  to  Xq/(K  - 1  )d,  the  angular  resolution  of  the  DPCA.  Both  CRLBs  are 
inversely  proportional  to  the  square  root  of  the  number  BW  of  independent  temporal 
samples  used  in  the  estimation  (for  statistically  homogeneous  reverberation)  and  to 
that  of  the  number  K  of  independent  elements  in  the  DPCA.  They  both  decrease 
with  the  reverberation-to-noise  ratio  p,  as  1/p  for  small  p  and  1  / yfp  for  large  p. 

The  above  expressions  of  the  CRLBs  show  the  benefits  of  operating  with  wide  band¬ 
width  B  for  DPCA  micronavigation,  as  both  DPCA  sway  and  yaw  accuracy  increase 
in  proportion  to  V BW  due  to  the  increase  in  the  number  of  independent  temporal 
samples.  In  addition  the  yaw  accuracy  is  seen  to  increase  in  proportion  to  A"3/2,  due 
to  the  combined  effect  of  the  increase  (as  K)  in  angular  resolution  of  the  DPCA  and 
the  increase  (as  Kl/ 2)  of  the  number  of  independent  spatial  samples.  Since 

K  -  U  7  2D~  (16) 

Lt 

K  can  be  increased  by  increasing  Lr ,  reducing  D  or  reducing  Lt,  i.e.  broadening  the 
transmission  beam.  Ultimately,  the  transmission  beam  will  be  limited  by  condition 
(5). 

The  bounds  (14)  and  (15)  were  first  presented  in  the  open  literature  by  Billon  and 
Fohanno[6]. 
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4.2  Experimental  validation 

It  was  shown  in  [8]  that  the  experimental  standard  deviations  approached  these 
CRJLBs  closely,  once  the  effect  of  calibration  errors  of  the  physical  array  were  re¬ 
moved.  A  more  complete  experimental  validation  is  presented  in  Fig.  4. 


x  io  3 


ping  pair 

Figure  4  Comparison  between  the  theoretical  CRLB  for  a  single  DPCA  element 
(red)  and  the  experimental  sample  standard  deviation  S  (black). 


The  data  were  collected  in  November  2000  at  Elba  Island  (Italy)  with  a  24  m  under¬ 
water  rail  deployed  4  m  above  a  seafloor  of  sand  and  rocks.  The  water  depth  varied 
between  13  m  and  15  m.  The  transmitter  of  length  Lt  =  2  cm  had  a  horizontal 
beamwidth  of  28  deg  and  a  vertical  beamwidth  of  2.8  deg.  The  transmitted  signal 
was  a  4  ms  LFM  swept  from  118  kHz  to  182  kHz,  with  a  ping  repetition  period 
of  0.5  s.  The  receive  array  had  32  elements  spaced  at  0.834  cm  for  a  total  length 
Lr  —  26.7  cm.  A  constant  along-track  velocity  of  6.675  cm/s  was  imposed  along  the 
rail,  giving  a  SAS  spatial  sampling  factor  a  =  4  and  24  elements  in  the  DPCA. 

The  DPCA  was  undcrsampled  by  a  factor  3,  to  arrive  at  8  independent  elements. 
Estimates  Jk  oi  n/k  {k  =  1, . . . ,  8)  were  obtained  using  a  short  term  correlation  with 
a  5  m  long  window  centred  at  45  m  range,  corresponding  to  W  —  6.7  x  10-3  s 
and  BW  =  427.  The  reverberation  to  noise  ratio  p,  computed  from  the  correlation 
coefficient  p  as  p  —  p/(l  -  p)  varies  between  4.9  dB  and  9.5  dB,  depending  on  the 
ping  pair  and  the  DPCA  element. 

In  Fig.  4  the  CRLB  (14)  is  plotted  as  a  function  of  the  ping  pair,  using  BW  =  427, 
K  =  1  and  the  p  averaged  over  the  8  DPCA  elements.  Under  the  assumption  of  no 
yaw  during  the  displacement  along  the  rail,  it  follows  from  Eq.  (3)  that  71, ...  ,7s 


-9- 


SACLANTCEN  SR-355 


should  be  identical  up  to  estimation  errors,  so  that  the  sample  standard  deviation  S 
provides  an  estimate  of  <j7.  S'  is  a  random  variable  which,  for  Gaussian  errors,  follows 
a  x-distribution  with  7  degrees  of  freedom.  It  is  shown  in  Fig.  4  that  S  is  in  good 
agreement  with  theory.  The  mean  value  of  the  relative  estimation  error  (S-a7)/tr7, 
obtained  averaging  over  the  ping  pairs,  is  only  about  6%.  Its  standard  deviation 
is  about  29%,  which  compares  well  with  the  theoretical  value  of  26%  derived  from 
the  x-dist.ribution.  The  accuracy  achieved  is  of  the  order  of  Ao/500,  which  is  quite 
remarkable. 
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5 

Accuracy  requirements  for  DPCA  sway  and  yaw 


Let  (0,x,y)  be  the  slant  range  plane,  with  Ox  along-track  and  Oy  across-track, 
(xp,yp)  the  coordinates  of  Cp  =  (Tp  +  Rp)/ 2,  where  Rp  is  the  centre  of  the  physical 
receive  aperture  at  ping  p,  9p  the  angle  between  Oy  and  boresight  to  the  physical 
aperture. 

The  trajectory  followed  by  the  sonar  can  be  expressed  as 


xp+ 1 

=  xp  +  D 

Up+ 1 

=  Pp  T  7p  T  ^2@p  T 

(17) 

Qp+ 1 

—  @p  + 

where  jp  and  typ  are  the  DPCA  sway  and  yaw  between  pings  p  and  p  +  1  and  the 
angles  9P  have  been  assumed  small  (Fig.  5). 

The  quantity  yp+\  -  yp,  which  can  be  interpreted  as  the  physical  sway  between 
successive  pings,  is  seen  to  be  the  sum  of  three  terms.  The  first  is  the  DPCA  sway 
and  the  other  two  result  from  the  heading  of  the  physical  reception  antenna  at  pings 
p  and  p+  1.  Indeed  the  geometrical  centre  of  the  DPCA  and  that  of  the  physical 
array  are  separated  by  D/2,  and  this  lever  arm  leads  to  a  difference  in  cross-track 
position  of  the  associated  phase  centres  of  D(9p  +  9p+ 1)/2  (Fig.  5).  After  algebraic 


Figure  5  SAS  trajectory  in  the  slant-range  plane. 
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manipulations,  the  estimated  trajectory  may  be  expressed  as: 
xp  =  ( p  -  1  )D 

<  yP  =  Tfc}'n  +  DZ&(p-i-i)ln  (is) 

1  0P  = 

It  is  important  that  this  trajectory  be  estimated  with  sufficient  precision.  Errors 
<$7p  and  5ipp  on  the  DPCA  estimates  will  accumulate  to 

IK  =  EClhi  +  DZti(p-‘-i)K  (19) 

K  =  miK 


The  most  important  effect  for  SAS  are  the  cross-track  micronavigation  errors  5yp. 
It  is  seen  in  the  first  equation  of  (19)  that  they  depend  on  accumulated  errors  on 
the  DPCA  sway  and  yaw.  In  the  presence  of  DPCA  sway  errors  alone  (Sipp  =  0) 
they  accumulate  like  a  random  walk,  whereas  in  the  presence  of  DPCA  yaw  errors 
alone  (Sjp  =  0)  they  accumulate  like  an  integrated  random  walk.  In  the  second  case 
the  errors  accumulate  much  faster  and  lead  to  a  highly  correlated  pattern  of  phase 
errors  along  the  SAS.  The  quantitative  analysis  of  these  errors  is  presented  in  what 
follows. 


5.1  Beam  pattern  specification  of  SAS  array  factor 


An  accuracy  study  was  presented  in  [6]  based  on  the  criterion  that  the  cross-track 
micronavigation  error  at  the  extremity  of  the  SAS,  6yp,  be  smaller  than  A/8.  This 
overlooks  the  fact  that  linear  phase  errors  steer  the  SAS  array  factor,  without  broad¬ 
ening  the  SAS  mainlobe,  so  that  the  term  in  the  8yp  which  is  linear  with  p  must 
be  removed  before  applying  the  criterion.  After  this,  the  criterion  should  be  similar 
to  the  standard  one  used  in  SAR  for  phase  errors  which  vary  quadratically  over 
the  aperture,  which  states  that  a  90  deg  phase  difference  between  the  centre  of  the 
aperture  and  its  extremities  leads  to  a  10%  increase  in  the  3  dB  beamwidth. 


However,  it  follows  from  (19)  that  if  the  Syp  were  expanded  as  function  of  p,  there 
would  certainly  be  terms  of  order  higher  than  2  which  are  not  accounted  for  in  [6]. 
A  more  appropriate  criterion  can  be  defined  using  the  normalized  gain  of  the  SAS 
array  factor  which  can  be  expressed  as 


9 


1 

P2 


Y1  exp(j2k0(8yp  -  ap )) 

p=l 


(20) 
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where  ko  —  27r/Ao  is  the  wavenumber  at  the  centre  frequency  and 


P 


°  =  £ 

p= i 


12(p-(P+l)/2) 
P(P2  -  1) 


(21) 


is  the  slope  of  the  linear  fit  of  8yp.  This  criterion  extends  that  of  [12]  where  the 
linear  term  had  been  overlooked. 


Subtracting  the  linear  fit  accounts  for  the  fact  that  cross-track  micronavigation 
errors  which  grow  linearly  with  p  simply  steer  the  SAS  array  factor  by  a/D,  without 
defocusing.  However,  the  steering  of  the  SAS  array  factor  away  from  the  pointing 
direction  of  the  physical  aperture  leads  to  an  increase  in  azimuth  ambiguities  as  well 
as  a  loss  in  reverberation  to  noise  ratio.  The  peak  of  the  grating  lobes  will  be  20  dB 
below  the  mainlobe  (see  [13]  eq.  (16)),  provided  that  a/D  <  Aq/20 D,  i.e., 


2a  <  l_ 
A0  ~  10' 


(22) 


The  SAS  beampattern  specification  will  be  given  in  terms  of 


9=  (9) 


(23) 


the  expected  normalized  gain  of  the  SAS  array  factor.  For  what  follows  it  is  also  con¬ 
venient  to  define  a  logarithmic  gain  G  =  10  log10  g  and  correspondingly  the  expected 
value  G.  Since  0  <  g  <  1,  one  has  G  <  0. 


By  elementary  algebraic  manipulation,  Eq.  (20)  becomes 

1  o  p  P_1 

9  =  -p  +  p2  Y,  £  cos(2 k0(8yp  -  Syq  -  a(p  -  q)))  (24) 

P=  1  9=1 


For  small  errors 
expected  gain 


(ko8yp  <$:  1)  Eq.  (24)  can  be  expanded  in  series  obtaining  the 


i  P  v- 1 


p2 


P=  1  9=1 


(25) 


where 

a2p<q  =  (\2ho{8yp  -  5yq  -  a(p  -  g))|2^  .  (26) 


Since  5jp  and  8tpp  are  independent,  the  corresponding  losses  in  SAS  array  gain  will 
be  studied  separately  in  the  next  two  subsections. 
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5.2  Required  DPCA  sway  accuracy 


We  assume  that  Sxjjp  —  0  for  all  p.  Then  by  (19)  one  has 

p- 1 

&Vp  -*Vq  =  Y,  (27) 

l=q 


and  substituting  (18)  into  (21) 

Q  =  p(p2'-T)^p(P~P)g> 


(28) 


After  algebraic  manipulations  of  (25),  one  finally  obtains 

<29) 

where  a 2  =  (Sj2).  The  required  accuracy  on  the  DPCA  sway  can  therefore  be 
expressed  as  a  function  of  P  and  g  as: 


cr7  < 


Ao  / 15(1  —  g) 
2  V  47 r2F 


(P»l). 


(30) 


It  is  seen  that,  for  given  g ,  the  standard  deviation  of  the  DPCA  sway  estimation  is 
required  to  decrease  as  P~1/2  for  large  P  (Fig.  6(a)). 


(b) 


Number  of  pings  P 


Figure  6  Accuracy  required  on  DPCA  sway  (a)  and  DPCA  yaw  (b),  to  sum  co¬ 
herently  P  pings  in  the  SAS,  for  given  values  of  the  normalized  expected  SAS  array 
gain  G. 
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Finally,  it  follows  from  (28)  that 

°a  =  ( P  »  !)•  (31) 

where  a2  =  (a2).  When  (30)  is  substituted  into  (31),  2cra/Ao  decreases  as  P-1  for 
large  P.  Thus,  condition  (30)  is  stronger  than  (22)  and  one  is  justified  in  specifying 
the  SAS  beampattern  solely  in  terms  of  the  SAS  array  factor. 


5.3  Required  DPCA  yaw  accuracy 


We  assume  Sjp  =  0  for  all  p.  Then  by  (19)  one  has 

Syp  ~6yq  =  D  J2(P  ~  q)Hi  +  D  |> -  l  -  \)5^ 

/=  1  l=q  1 


and 


n^(P2-l+pP-2p2)(P-p) 
a  =  D2^ - mP2_n - Wp 

P=  1 


P(P2  -  1) 

After  rather  cumbersome  algebraic  manipulations,  g  can  be  expressed  as 


(32) 


(33) 


(34) 


where  —  (Sift2). 

The  required  accuracy  for  the  yaw  can  therefore  be  expressed  as  a  function  of  P  and 
9  as: 


Ap  Im(l-g) 

—  2D  V  7 r2P3  (”»!)•  (33) 

The  standard  deviation  of  the  DPCA  yaw  estimation  is  now  required  to  decrease 
much  faster  with  P  than  for  the  DPCA  sway,  in  proportion  to  P-3/2  for  large  P 
instead  of  P-1/2  (Fig.  6(b)). 

Finally,  it  follows  from  (33)  that 


aa  =  \j  (P  >  (36) 

When  (35)  is  substituted  into  (36),  2<7a/Ao  decreases  as  P-1  for  large  P.  Thus,  as 
before,  condition  (35)  is  stronger  than  (22). 
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6 

Required  reverberation-to-noise  ratio 


By  combining  the  results  obtained  in  Sections  4  and  5  it  is  now  possible  to  find 
the  sonar  performance  requirements  needed  to  attain  a  given  Q.  To  facilitate  the 
interpretation,  it  is  useful  to  introduce  the  effective  reverberation  to  noise  ratio 


1  ■+■  2  p 


(37) 


For  p  >  1  this  can  be  interpreted  as  the  product  of  the  reverberation-to-noise  ratio 
p,  the  number  Lr/Lt  of  independent  elements  in  the  receive  aperture  and  the  number 
BW  of  independent  temporal  samples. 

The  CRLBs  (14)  and  (15)  of  the  DPCA  sway  and  yaw  estimators  can  be  then 
expressed  respectively  by 


1  Ao 

1  « 

1 

(38) 

27 r  2  1 

J  a-l  y/p^ 

y/3  Ao 

r « 

1 

(39) 

7T  2D  \l 

'  (a  -  1)3 

where  in  the  latter  it  has  been  assumed  K2  1. 


By  substituting  (38)  and  (39)  into  (30)  and  (35)  it  is  straightforward  to  obtain  the 
effective  reverberation  to  noise  ratios  required  to  limit,  respectively,  the  DPCA  sway 
and  yaw  errors  as  a  function  of  Q,  a  and  1  —  g.  In  the  limit  P  1,  for  which  (7) 
becomes  Q  -  1  =  P/a ,  they  read 


Pea  P 


Q  -  1  a2 
15(1  —  g)  a  —  1 


(40) 


>  (<?  ~  l)3  «4 

P°«-  35(1-5)  (<* -l)3 


(41) 


where  it  has  been  assumed  that  1  —  g  <C  1.  The  required  pefi  established  in  [12]  was 
about  11.5  dB  higher,  due  to  the  omitted  linear  term  in  (20). 
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(a)  requirements  for  sway  (b)  requirements  for  yaw 


Figure  7  Required  peff  for  DPCA  sway  (a)  and  DPCA  yaw  (b),  estimated  as  a 
function  of  desired  resolution  gain  Q  and  SAS  sampling  factor  a,  for  a  normalized 
expected  SAS  array  gain  G  =  —0.25  dB. 


In  Fig.  7  the  values  of  peff  needed  to  obtain  G  >  —0.25  dB  are  plotted  as  functions  of 
Q  and  a  for  the  two  cases.  Figure  7(a)  represents  the  requirements  in  the  limit  when 
only  the  DPCA  sway  has  to  be  estimated  (i.e.  for  an  infinitely  precise  independent 
estimate  of  the  yaw)  whereas  Fig.  7(b)  represents  those  requirements  when  both 
DPCA  sway  and  yaw  have  to  be  estimated. 

The  requirements  for  accurate  DPCA  yaw  estimation  are  seen  to  be  higher,  by 
several  orders  of  magnitude,  than  those  for  the  DPCA  sway.  Thus  the  achievable  Q 
by  a  DPCA  micronavigated  SAS  is  limited  by  estimation  errors  on  the  heading  of 
the  physical  array,  which  induce  cross-track  position  errors  of  the  physical  array  in 
DPCA  micronavigation. 

The  right  hand  side  of  (41)  attains  an  absolute  minimum  for  a  =  4.  Therefore,  a 
tradeoff  between  resolution  gain  and  area  mapping  rate  characterizes  the  design  of  a 
DPCA-micronavigated  SAS.  According  to  (41)  this  tradeoff  will  usually  be  in  favor 
of  a  <  4.  For  fixed  peff,  the  Q  obtained  with  a  =  2  has  a  loss  of  only  16%  with 
respect  to  the  optimal  a.  The  loss  is  much  more  severe  for  a  approaching  1. 


6.1  Simulation 


To  illustrate  the  preceding  analytical  results,  various  realizations  of  the  5yp  and 
corresponding  beampatterns  were  simulated  for  Q  =  10,  a  =  4  and  peff  =  35.4  dB, 
which  correspond  to  a  gain  G  =  —0.25  dB  according  to  (41).  The  first  five  realiza- 
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tions  are  illustrated  in  Fig.  8,  followed  by  the  reference  case  with  no  micronavigation 
errors.  The  8yp  were  obtained  by  generating  independent  normal  errors  Sipp  of  zero 
mean  and  standard  deviation  given  by  the  approximate  CRLBs  (39)  for  the  above 
values  of  a  and  peff  and  summing  them  according  to  the  first  equation  in  (19).  It 
should  be  noted  that  the  8yp  have,  as  expected,  a  high  degree  of  correlation. 


Figure  8  Simulated  cross-track  micronavigation  errors  8yp  (a),  corresponding  SAS 
array  factor  beampatterns  (b)  and  full  SAS  beampatterns  (c),  for  Q  =  10,  a  —  4  and 
ptff  —  35.4  dB.  In  (b)  are  also  indicated  the  normalized  SAS  array  factor  gains 
obtained  in  the  realizations. 


The  expected  SAS  gain,  estimated  by  means  of  (20)  by  averaging  over  106  realiza¬ 
tions,  is  equal  to  G  =  —0.235  dB,  in  close  agreement  with  the  theoretical  value  of 
G  =  —0.25  dB.  The  small  discrepancy  is  attributed  to  approximation  (25).  The 
steering  of  the  SAS  array  factor  is  apparent  in  Fig.  8(b),  but  it  is  seen  to  be  of  no 
consequence  for  the  full  SAS  beampattern,  as  no  grating  lobes  are  visible  at  ±1  in 
Fig.  8(c). 
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The  simulation  also  allows  the  estimation  of  the  whole  distribution  of  G  for  given 
Q ,  a  and  peff.  For  Q  =  10,  a  =  4  and  peff  =  35.4  dB,  the  gain  G  was  found  larger 
than  -0.84  dB  for  95%  of  the  realizations. 


6.2  Comparison  with  experiments 

DPCA-micronavigated  SAS  results  for  a  dual  frequency  405/165  kHz  sonar  were 
presented  in  [6].  The  experimental  parameters  of  the  405  kHz  SAS  were  a  =  1.7, 
BW  =  156,  Lr/Lt  =  29,  p  =  7.5  dB.  One  has  then  peff  —  44  dB  so  that  (41),  with 
G  =  —0.25  dB,  gives  a  theoretical  Q  of  14.  According  to  [6],  the  Q  achieved  experi¬ 
mentally  was  close  to  9,  although  a  direct  experimental  measure  of  the  beampattern 
was  not  possible  due  to  the  absence  of  a  strong  point  scatterer. 

For  the  165  kHz  SAS,  the  experimental  parameters  were  a  =  1.8,  BW  —  156, 
Lr/Lt  =  20,  p  =  6  dB.  One  has  then  pen  ~  40.9  dB  giving  now  a  theoretical  Q 
of  11.5.  According  to  [6],  the  Q  achieved  experimentally  is  close  to  12.  However, 
the  experimental  beampattern,  measured  for  a  spherical  target,  shows  peak  sidelobe 
levels  only  4  dB  below  the  main  peak. 

For  the  100  kHz  DPCA-micronavigated  SAS  discussed  in  [9],  the  data  relative  to 
Fig.  6  were  obtained  with  BW  =  266,  Lr/Lt  =  45.6,  p  =  9.5  dB  and  a  =  4.3,  giving 
peff  =  50.4  dB  and  a  theoretical  Q  of  29  to  be  compared  to  the  value  of  Q  =  11 
obtained  experimentally. 

The  above  results  show  that  the  SAS  performance  obtained  experimentally  is  usu¬ 
ally  significantly  lower  than  the  optimum  theoretical  performance,  determined  by 
the  CRLBs.  The  origins  of  the  discrepancy  are  as  yet  unexplained.  Amongst  the 
possible  causes  are  residual  calibration  errors  of  the  physical  elements  as  well  as 
approximations  such  as  (5)  and  to  a  lesser  extent  (2). 

The  design  of  a  180/20  kHz  DPCA-micronavigated  SAS  is  presented  in  [5].  For  the 
180  kHz  HFSAS  the  receive  array  is  comprised  of  11  elements  spaced  at  5  cm  and 
the  desired  cross-range  resolution  up  to  a  maximum  range  of  40  m  is  2.5  cm,  i.e., 
Q  =  24  at  the  far  range.  For  the  20  kHz  LFSAS  the  receive  array  is  comprised  of  14 
elements  spaced  at  7.5  cm  and  the  desired  cross-range  resolution  upto  40  m  is  7.5 
cm,  i.e.  Q  =  76  at  the  far  range.  Both  SAS  are  designed  with  only  two  elements  in 
the  DPCA,  i.e.,  K  =  2,  giving  a  —  14/12  for  the  180  kHz  HFSAS  and  a  —  11/9  for 
the  20  kHz  LFSAS.  It  follows  then  from  (41)  that  the  required  peff  is  of  the  order  of 
64.4  dB  for  the  HFSAS  and  76.5  dB  for  the  LFSAS. 

Clearly  obtaining  such  high  values  is  challenging.  The  good  experimental  results 
obtained  may  be  due  to  the  fact  that  the  HF/LFSAS  is  mounted  on  a  very  stable 
towbody.  Assuming  the  towbody  keeps  a  constant  heading  during  the  SAS  integra- 
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tion  time,  the  assumption  that  there  is  no  yaw  is  valid  so  that  (40),  rather  that  (41), 
applies.  The  requirements  in  terms  of  peS  are  then  much  more  modest,  of  the  order 
of  20.3  dB  for  the  HFSAS  and  20.7  dB  for  the  LFSAS. 
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7 

Conclusions 


The  two  most  important  parameters  which  govern  the  cross-range  resolution  gain  Q 
achievable  by  a  DPCA-micronavigated  SAS  are: 


1.  the  spatial  oversampling  factor  a, 

2.  the  effective  reverberation  to  noise  ratio  peff,  defined  as  the  product  of  the 
reverberation-to-noise  ratio  p,  the  number  Lr/Lt  of  independent  elements  in 
the  receive  aperture  and  the  number  BW  of  independent  temporal  samples 
used  in  the  estimation. 


The  achievable  Q ,  for  given  peff,  can  be  maximized  by  operating  close  to  the  optimum 
value  of  a  which  was  shown  to  be  a  =  4.  The  price  to  pay,  however,  is  a  reduction 
by  a  factor  a  of  the  area  mapping  rate. 

The  achievable  Q  is  limited  chiefly  by  estimation  errors  on  the  heading  of  the  physical 
array,  which  induce  cross-track  DPCA  micronavigation  errors  of  the  SAS.  Means 
of  aiding  DPCA  in  estimating  this  heading,  such  as  additional  instrumentation  or 
new  algorithms,  would  make  possible  higher  resolution  gains  and  increased  mapping 
rates. 
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